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Transient and Steady-State Velocity of Domain Walls 
for a Complete Range of Drive Fields 

Henry C. Bourne, Jr., and David S. Bartran 

Abstract - Approximate analytic solutions for transient and steady-state 
180® domain wall motion in bulk magnetic material are ^obtained from the 
dynamic torque equations with a Gilbert damping term. The results for the 
Walker region in which the transient solution approaches the familiar Walker 
steady-state solution are presented in a slightly new form for completeness. 

An analytic solution Corresponding to larger drive fields predicts an oscil- 
latory motion with an average value which decreases with drive field for 
reasonable values of the damping parameter. These results agree with those 
obtained by a computer solution of the torque equation and those obtained with 
the assumption of a very large anisotropy field. 


INTRODUCTION 


( 1 ) 


In a previous paper'* , transient solutions for the motion. of domain 
walls in bulk magnetic materials were obtained from the vector equation of 
motion, 

— = X (H - — — a) 
in which M is the saturation magnetization, y is the gyromagnetic ratio, 
H is the effective field including applied, stray, anisotropy, and exchange 
components, and a is the damping parameter. With the coordinate system 
and wall configuration as represented in Fig. 1, the solutions obtained for 
e, q), and the corresponding wall velocity, v, are as follows 


This work was supported by NSF Grant GH 34584 and NASA Grant NGR-44-006- 
001. H. Bourne is with the Department of Electrical Engineering, Rice Univ- 
ersity, Houston, Texas. D. Bartran is with Honeywell, Oklahoma City, Oklahoma. 
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w 

'C.n(tan 0 / 2 ) » C 2 <t)[y - J v(z)dz} 

0 


< 2 ) 


^ Vis ‘p + V 
0 


and 


V » Y-,/"^ 

p, 


1/2 


orh^ + sin tp cos cp 

o (I4<y^)(8in^<p+hj^>^^^ 


in which cp is obtained from 


(3) 


(4) 


2 h 


(l + a^) =f ^ - sin 2 9 ) . 


(5) 


In these equations A is the exchange constant, h^^ is the normalized anisotropy 

-7 

field p, H- /M , h is the normalized applied field u H /M, and a = 4 tt X 10 mks 
o k ' z *^o z * ^o 

The only approximation Involved concerns a consistency condition which requires 


that 

— [y - J v(T)dT] « C^(t) v(t) 

in order for cp ^ f (y) , a basic assumption in the trial solution. The solution 
under this assumption indicates that this approximation is an excellent one 
for the wall widths encountered in magnetic materials. This condition is 
equivalent to neglecting the incremental velocity along the wall compared to 
the velocity of the wall center. An alternate assumption which gives the same 
result concerns the partition of Eq. 1 such that the same terras equate to 
determine the structure constant, C^, and the velocity at the wall center, v, 
during the transient as in the steady^state solution, which involves no ap- 
proximation at least In the region 2 h /<y < 1 , 

2 

For a stationary wall h 9 « 0 and the equations reduce to the familiar 

z 


3 


'tn(tan 9/2) = y '^ k7a , For 2h^/or < 1, the Walker solution^^^ corresponding 

to 2 h^/o' « sin 2 cp predicts a contracted wall moving at constant velocity. 

The transient solution corresponding to this constant velocity solution is 

the one examined theoretically and experimentally in previous papers 

For completeness and for comparison, this solution is repeated in a slightly 

different form together with the solution for 2h /oi > 1. All solutions 

2 

assume a step function drive field, h . 

2 

Walker Region; 2h^/of < 1 

Eq. 5, which predicts that <p increases monotonlcally In time until it 
reaches a steady-state value corresponding to 2 h /or * sin 29 , is integrated 
to yield 

2h^/or 

' tan 9 « 

with 

9=^1 •\j\ - (,2-hJa)^ /(I + 

If the appropriate trigonometric functions of 9 obtained from Eq, 6 are 

substituted in Eq. 4 and the hyperbolic function is expressed in terms of 

exponentials, the ’’exact’* solution of reference (1) is obtained. Note 

that a factor in the effective damping 0 Involves the drive field. With 

Ihjoi < 1 , the steady-state value of 9 is less than tt/ 4. Although 9 itself 

is well-behaved, Eq. 4 is sufficiently complex that the transient response 

contains an overshoot, for 2 h,^^^(l -h h )^^^f(l + h. ^ < 1 

k tc k k Of • * 

as previously shown. 

Limiting Case: 2h /of « 1 

2 


1 +*yi-(2h^/a) 


i 
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2h 

In the limiting case, " I, either from Eq. 6 or directly from 

integration of Eq. 5, 

(\M/p, )oitf2 

tan cp “ ■ , , 1 .... — 

(YM/p,^)0ft/2 + (1 + Of^) 

The angle cp increases monotonically to n/4 and the solution is similar to 

that of the previous case, with an overshoot in the transient response. 

Oscillatory Region; 2h /(V > 1 

2 

In this particular case Eq, 5 predicts that q> continues to increase for 

all time with a periodically changing rate. From Eq. 4 the velocity behaves 

as sin 2cp displaced by the relatively small terra ofh and modified by a posi- 

z 

tive term of oscillating magnitude in the denominator so that the wall velocity 
is alternately positive (0 < q> < tt/2) and negative (tt/ 2 < q> < tt) and repeats. 
Again Eq. 5 may be Integrated directly to obtain 

2h foi 
z 

tan <p.s» ■ ' — 

- 1 cot U)t 

with 

(s = - ^(2h^/a)^ -1 /(I + Qf^) 

The periodic nature of the solution is apparent. The system is in a 

steady-state oscillatory condition from the beginning with no transient in« 

volved. The frequency depends on the drive field, h . In the limit of 

z 

^h /o = 1 the previous solution is obtained. 

2 * 

The maximum velocity magnitude in time is the same as the maximum steady- 
state velocity in the Walker region and corresponds to 


5 


tan 9 « ± /(I + , 

sin cp cos (j) =,+ + hj^)’-^^((l + , 

sin% = , 


and 


max 


= vi/5 ta + h^)^/2 - , 


independent of the drive field, h . However, the time at which the maximum 

2 

occurs depends on h and h and is given by 

» z 


cot U)t 


- 1 +(2h^/a)(l + 
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Figs. 2, 3, and 4 give the velocity as a function of time for several values 

of drive field and three values of h^^ corresponding to permalloys and 

bubble-type materials. 

2 21^2 

For or (r^) « the velocity is approximately equal to zero for 
9 = 0,TT which corresponds to tnt « 0, tt , The velocity is also approximately 
equal to zero for cp = tt/ 2 which corresponds to 


cot U)t I 

v=0 


1 

■y(2h^/Q')2 - 1 


In general more time is spent in the positive velocity region and the velocity 
has a nonzero average value which is shown in Fig. 5. 

j 2h 

In the limit of very large drives, 2h^/of » 1 , but again « 1, 

the positive and negative half-cycles of velocty are an odd function about 

oJt = tt/ 2, cp ^ u)t Y H t and the average velocity is zero. If in addition 

z 

hj^ » 1 the velocity behaves as sin 2u>t at very high frequencies. 
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Conclusioa 

With suitable approximations a complete solution of the vector equation 
tof motion with a viscous damping parameter may be obtained which describes 
the motion of a magnetic domain wall over a comple« range of drive fields. 
The solution pertains to a bulk magnetic material 'scribed by an exchange 
constant, an anisotropy constant, a saturation magnetization, and a damping 
parameter. The oscillatory response predicted in the high drive region 
probably may best be explored experimentally in materials with relatively 


small saturation magnetizations and large anisotropy constants. 

These results seem to agree with those reported by Slonczewski^^^ and 
/ 6 \ 

by Walker and Schryer^ in a computer solution of the torque equation. 
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Figure Captions 


Fig. 1. Coordinate System and Wall Configuration 

Fig. 2. Normalized Velocity as a Function of Time for h, « u H, /M « 0.0005 

k ok 

and Various Drive Fields 2h /of = 2u, H /orM 

z z 


Fig, 3. Normalized Velocity as a Function of Time for h. 

k 

and Various Drive Fields 2h /a » 2u- H '/oM 

z ^o z 

Fig. 4. Normalized Velocity as a Function of Time for h. 

and Various Drive Fields 2h /(y = 2u H /aM 

2 *^0 2 
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Fig. 5. Average Velocity as a Function of Drive Field for Various h. 

k 

in the Walker Region and in the Oscillatory Region 
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